Integrative and Comparative Analysis of Retinoblastoma and Osteosarcoma

In the last one and a half decades, the generalization of high throughput methods in molecular biology has led to the generation of vast amounts of datasets that unraveled the unfathomed complexity of the cell regulatory mechanisms. The recently published results of the ENCODE project (ENCODE Project Consortium et al., 2012) demonstrated the extend of these in the human genome and certainly more regulation mechanisms will be discovered in the future. Already, this complexity within a single cell - without taking into account cell-cell interaction or micro-environment influences - cannot be abstracted by the human mind. However, understanding it is the key to devise adapted treatments to genetic diseases or disorders, among which is cancer. In mathematics, such complex problems are addressed using methods that reduce their complexity, so that they can be modeled in a solvable manner. In biology, it led researchers to develop the concept of systems biology as a mean to abstract the complexity of the cell regulatory network. To date, most of the published studies using high throughput technologies only focus on one kind of regulatory mechanism and hence cannot be used as such to investigate the interactions between these. Moreover, distinguishing causative from confounding factors within such studies is difficult. These were my original motivations to develop analytical and statistical methods that control for confounding factors effects and allow the integrative and comparative analysis of different kinds of datasets. In fine, three different tools were developed to achieve this goal. First, "customCDF": a tool to redefine the Custom Definition File (CDF) of Affymetrix GeneChips. It results in the increased sensitivity of downstream analyses as these bene fit from the constantly evolving human genome reference and annotations. Second, "aSim": a tool to simulate microarray data, which was required to benchmark the developed algorithms. Third, for the integrative analysis, a set of combined statistical methods and finally for the comparative analysis, a modification of the integrative analysis approach. These were bundled in the "crossChip" R package. The "customCDF" and "aSim" tools were first validated on independant datasets. The developed analytical methods ("crossChip") were first validated on "aSim" simulated data and publicly available datasets and then used to answer two biological questions. First, using two retinoblastoma datasets, the effect of genomic copy number variations on gene-expression was investigated. Then, motivated by the fact that retinoblastoma patients have a higher chance to develop osteosarcoma later in life than the average population, datasets of both these tumors were comparatively analyzed to assess these tumors similarities and differences. Despite a rather limited number of samples within the selected datasets, the developed approaches with their higher sensitivity and sensibility were successful and set the ground for larger scale analyses. Indeed, the integrative analysis applied to retinoblastoma revealed the high importance of the chromosome 6 gain at a later stage of the disease, indicating that many genes on that chromosome are beneficial to cancerogenesis. Moreover, in comparison to standard microarray analyses, it demonstrated its efficacy at detecting the interplay of regulatory mechanisms: examples of positive and negative compensation of gene expression in lost and gained regions, respectively, as well as examples of antisense transcription, pseudogene and snRNAs regulation were identified in this dataset. The comparative analysis on the other hand revealed the high similarity of the retinoblastoma and osteosarcoma tumors, while at the same time showing that either of them take advantage of their distinct micro-environment and consequently appear to make use of different signaling pathways, PKC/calmodulin in retinoblastoma and GPCR/RAS in osteosarcoma. The developed tools and statistical methods have demonstrated their validity and utility by giving sensible answers to the two biological questions addressed. Moreover, they generated a large number of interesting hypotheses that need further investigations. And as they are not limited to microarray analysis but can be applied to analyze any high-throughput generated data, they demonstrated the usefulness of "systems biology" approaches to study cancerogenesis

Exploring the beta-tubulin gene family in a benzimidazole-resistant Parascaris univalens population

Benzimidazole (BZ) drugs are frequently used to treat infections with the equine ascarid Parascaris univalens due to increasing resistance to macrocyclic lactones and pyrantel. Benzimidazole resistance is rare in ascarids in contrast to strongyle parasites where this resistance is widespread. In strongyles, single nucleotide polymorphisms (SNPs) at codons 167, 198 and 200 in a 13-tubulin gene have been correlated to BZ resistance, but little is known about the 13-tubulin genes and their possible involvement in BZ resistance in P. univalens and other ascarids. Previously two 13-tubulin genes have been identified in P. univalens. In this study, we present five additional 13-tubulin genes as well as the phylogenetic relationship of all seven genes to 13-tubulins of other clade III and V nematodes. In addition, the efficacy of fenbendazole for treatment of P. univalens on a Swedish stud farm was studied in 2019 and 2020 using faecal egg count reduction test. Reductions varied from 73% to 88%, indicating the presence of a resistant P. univalens population on the farm. The emergence of BZ resistance emphasizes the need for development of molecular markers for rapid and more sensitive detection of resistant populations. We therefore investigated whether possible SNPs at positions 167, 198 or 200 in any of the 13-tubulin genes could be used to distinguish between resistant and susceptible P. univalens populations. Amplicon sequencing covering the mutation sites 167, 198 and 200 in all seven 13-tubulin genes revealed an absence of SNPs in both resistant and susceptible populations, suggesting that the mechanism behind BZ resistance in ascarids is different from that in strongyle nematodes and the search for a molecular marker for BZ resistance in P. univalens needs to continue

Training in High-Throughput Sequencing: Common Guidelines to Enable Material Sharing, Dissemination, and Reusability.

This is the author accepted manuscript. It is currently under an indefinite embargo pending publication by the Public Library of Science.The advancement of high-throughput sequencing (HTS) technologies and the rapid development of numerous analysis algorithms and pipelines in this field has resulted in an unprecedentedly high demand for training scientists in HTS data analysis. Embarking on developing new training materials is challenging for many reasons. Trainers often do not have prior experience in preparing or delivering such materials and struggle to keep them up to date. A repository of curated HTS training materials would support trainers in materials preparation, reduce the duplication of effort by increasing the usage of existing materials, and allow for the sharing of teaching experience among the HTS trainers' community. To achieve this, we have developed a strategy for materials' curation and dissemination. Standards for describing training materials have been proposed and applied to the curation of existing materials. A Git repository has been set up for sharing annotated materials that can now be reused, modified, or incorporated into new courses. This repository uses Git; hence, it is decentralized and self-managed by the community and can be forked/built-upon by all users. The repository is accessible at http://bioinformatics.upsc.se/htmr.COST action SeqAhead (BM1006

ProkSeq for complete analysis of RNA-Seq data from prokaryotes

Since its introduction, RNA-Seq technology has been used extensively in studies of pathogenic bacteria to identify and quantify differences in gene expression across multiple samples from bacteria exposed to different conditions. With some exceptions, tools for studying gene expression, determination of differential gene expression, downstream pathway analysis and normalization of data collected in extreme biological conditions is still lacking. Here, we describe ProkSeq, a user-friendly, fully automated RNA-Seq data analysis pipeline designed for prokaryotes. ProkSeq provides a wide variety of options for analysing differential expression, normalizing expression data and visualizing data and results

FACT – a framework for the functional interpretation of high-throughput experiments

BACKGROUND: Interpreting the results of high-throughput experiments, such as those obtained from DNA-microarrays, is an often time-consuming task due to the high number of data-points that need to be analyzed in parallel. It is usually a matter of extensive testing and unknown beforehand, which of the possible approaches for the functional analysis will be the most informative RESULTS: To address this problem, we have developed the Flexible Annotation and Correlation Tool (FACT). FACT allows for detection of important patterns in large data sets by simplifying the integration of heterogeneous data sources and the subsequent application of different algorithms for statistical evaluation or visualization of the annotated data. The system is constantly extended to include additional annotation data and comparison methods. CONCLUSION: FACT serves as a highly flexible framework for the explorative analysis of large genomic and proteomic result sets. The program can be used online; open source code and supplementary information are available at

Transcriptomics of ivermectin response in Caenorhabditis elegans: Integrating abamectin quantitative trait loci and comparison to the Ivermectin-exposed DA1316 strain

Parasitic nematodes pose a significant threat to human and animal health, as well as cause economic losses in the agricultural sector. The use of anthelmintic drugs, such as Ivermectin (IVM), to control these parasites has led to widespread drug resistance. Identifying genetic markers of resistance in parasitic nematodes can be challenging, but the free-living nematode Caenorhabditis elegans provides a suitable model. In this study, we aimed to analyze the transcriptomes of adult C. elegans worms of the N2 strain exposed to the anthelmintic drug Ivermectin (IVM), and compare them to those of the resistant strain DA1316 and the recently identified Abamectin Quantitative Trait Loci (QTL) on chromosome V. We exposed pools of 300 adult N2 worms to IVM (10(-7) and 10(-8) M) for 4 hours at 20 degrees C, extracted total RNA and sequenced it on the Illumina NovaSeq6000 platform. Differentially expressed genes (DEGs) were determined using an in-house pipeline. The DEGs were compared to genes from a previous microarray study on IVM-resistant C. elegans and Abamectin-QTL. Our results revealed 615 DEGs (183 up-regulated and 432 down-regulated genes) from diverse gene families in the N2 C. elegans strain. Of these DEGs, 31 overlapped with genes from IVM-exposed adult worms of the DA1316 strain. We identified 19 genes, including the folate transporter (folt-2) and the transmembrane transporter (T22F3.11), which exhibited an opposite expression in N2 and the DA1316 strain and were deemed potential candidates. Additionally, we compiled a list of potential candidates for further research including T-type calcium channel (cca-1), potassium chloride cotransporter (kcc-2), as well as other genes such as glutamate-gated channel (glc-1) that mapped to the Abamectin-QTL

Transcriptome Analysis of an Aedes albopictus Cell Line Single- and Dual-Infected with Lammi Virus and WNV

Understanding the flavivirus infection process in mosquito hosts is important and fundamental in the search for novel control strategies that target the mosquitoes' ability to carry and transmit pathogenic arboviruses. A group of viruses known as insect-specific viruses (ISVs) has been shown to interfere with the infection and replication of a secondary arbovirus infection in mosquitoes and mosquito-derived cell lines. However, the molecular mechanisms behind this interference are unknown. Therefore, in the present study, we infected the Aedes albopictus cell line U4.4 with either the West Nile virus (WNV), the insect-specific Lammi virus (LamV) or an infection scheme whereby cells were pre-infected with LamV 24 h prior to WNV challenge. The qPCR analysis showed that the dual-infected U4.4 cells had a reduced number of WNV RNA copies compared to WNV-only infected cells. The transcriptome profiles of the different infection groups showed a variety of genes with altered expression. WNV-infected cells had an up-regulation of a broad range of immune-related genes, while in LamV-infected cells, many genes related to stress, such as different heat-shock proteins, were up-regulated. The transcriptome profile of the dual-infected cells was a mix of up- and down-regulated genes triggered by both viruses. Furthermore, we observed an up-regulation of signal peptidase complex (SPC) proteins in all infection groups. These SPC proteins have shown importance for flavivirus assembly and secretion and could be potential targets for gene modification in strategies for the interruption of flavivirus transmission by mosquitoes

FLOWERING LOCUS T paralogs control the annual growth cycle in Populus trees

In temperate and boreal regions, perennials adapt their annual growth cycle to the change of seasons. These adaptations ensure survival in harsh environmental conditions, allowing growth at different latitudes and altitudes, and are therefore tightly regulated. Populus tree species cease growth and form terminal buds in autumn when photoperiod falls below a certain threshold.(1) This is followed by establishment of dormancy and cold hardiness over the winter. At the center of the photoperiodic pathway in Populus is the gene FLOWERING LOCUS T2 (FT2), which is expressed during summer and harbors significant SNPs in its locus associated with timing of bud set.(1-4) The paralogous gene FT1, on the other hand, is hyper-induced in chilling buds during winter.(3,5) Even though its function is so far unknown, it has been suggested to be involved in the regulation of flowering and the release of winter dormancy(.3,5). In this study, we employ CRISPR-Cas9-mediated gene editing to individually study the function of the FT-like genes in Populus trees. We show that while FT2 is required for vegetative growth during spring and summer and regulates the entry into dormancy, expression of FT1 is absolutely required for bud flush in spring. Gene expression profiling suggests that this function of FT1 is linked to the release of winter dormancy rather than to the regulation of bud flush per se. These data show how FT duplication and sub-functionalization have allowed Populus trees to regulate two completely different and major developmental control points during the yearly growth cycle

Comparative Fungal Community Analyses Using Metatranscriptomics and Internal Transcribed Spacer Amplicon Sequencing from Norway Spruce

The health, growth, and fitness of boreal forest trees are impacted and improved by their associated microbiomes. Microbial gene expression and functional activity can be assayed with RNA sequencing (RNA-Seq) data from host samples. In contrast, phylogenetic marker gene amplicon sequencing data are used to assess taxonomic composition and community structure of the microbiome. Few studies have considered how much of this structural and taxonomic information is included in transcriptomic data from matched samples. Here, we described fungal communities using both host-derived RNA-Seq and fungal ITS1 DNA amplicon sequencing to compare the outcomes between the methods. We used a panel of root and needle samples from the coniferous tree species Picea abies (Norway spruce) growing in untreated (nutrient-deficient) and nutrient-enriched plots at the Flakaliden forest research site in boreal northern Sweden. We show that the relationship between samples and alpha and beta diversity indicated by the fungal transcriptome is in agreement with that generated by the ITS data, while also identifying a lack of taxonomic overlap due to limitations imposed by current database coverage. Furthermore, we demonstrate how metatranscriptomics data additionally provide biologically informative functional insights. At the community level, there were changes in starch and sucrose metabolism, biosynthesis of amino acids, and pentose and glucuronate interconversions, while processing of organic macromolecules, including aromatic and heterocyclic compounds, was enriched in transcripts assigned to the genus Cortinarius.IMPORTANCE A deeper understanding of microbial communities associated with plants is revealing their importance for plant health and productivity. RNA extracted from plant field samples represents the host and other organisms present. Typically, gene expression studies focus on the plant component or, in a limited number of studies, expression in one or more associated organisms. However, metatranscriptomic data are rarely used for taxonomic profiling, which is currently performed using amplicon approaches. We created an assembly-based, reproducible, and hardware-agnostic workflow to taxonomically and functionally annotate fungal RNA-Seq data obtained from Norway spruce roots, which we compared to matching ITS amplicon sequencing data. While we identified some limitations and caveats, we show that functional, taxonomic, and compositional insights can all be obtained from RNA-Seq data. These findings highlight the potential of metatranscriptomics to advance our understanding of interaction, response, and effect between host plants and their associated microbial communities

PopulusPtERF85 Balances Xylem Cell Expansion and Secondary Cell Wall Formation in Hybrid Aspen

Secondary growth relies on precise and specialized transcriptional networks that determine cell division, differentiation, and maturation of xylem cells. We identified a novel role for the ethylene-induced Populus Ethylene Response Factor PtERF85 (Potri.015G023200) in balancing xylem cell expansion and secondary cell wall (SCW) formation in hybrid aspen (Populus tremula x tremuloides). Expression of PtERF85 is high in phloem and cambium cells and during the expansion of xylem cells, while it is low in maturing xylem tissue. Extending PtERF85 expression into SCW forming zones of woody tissues through ectopic expression reduced wood density and SCW thickness of xylem fibers but increased fiber diameter. Xylem transcriptomes from the transgenic trees revealed transcriptional induction of genes involved in cell expansion, translation, and growth. The expression of genes associated with plant vascular development and the biosynthesis of SCW chemical components such as xylan and lignin, was down-regulated in the transgenic trees. Our results suggest that PtERF85 activates genes related to xylem cell expansion, while preventing transcriptional activation of genes related to SCW formation. The importance of precise spatial expression of PtERF85 during wood development together with the observed phenotypes in response to ectopic PtERF85 expression suggests that PtERF85 contributes to the transition of fiber cells from elongation to secondary cell wall deposition